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Abstract. This manuscript is a survey on results related to boundary layers 
and the vanishing viscosity limit for incompressible flow. It is the lecture 
notes for a 10 hour minicourse given at the Morningside Center, Academia 
Sinica, Beijing, PRC from 11/28 to 12/07, 2007. The main topics covered are: 
a derivation of Prandtl's boundary layer equation; an outline of the rigorous 
theory of Prandtl's equation, without proofs; Kato's criterion for the vanishing 
viscosity limit; the vanishing viscosity limit with Navier friction condition; 
rigorous boundary layer theory for the Navier friction condition and boundary 
layers for flows in a rotating cylinder. 



1. Introduction 

In 1904, the issue of heavier-than-air, self propelled flight by human-made ma- 
chines was at the very edge of both science and technology. The first such flight, 
by the Wright brothers, occurred at December 14th, 1903. A Brazilian author is 
honor bound to remark that a more satisfying, and better publicized "first flight" 
was achieved by Santos Dumont, a Brazilian living in France, in September, 1906. 
The flight of a fixed-wing airplane could, at least in principle, be described by near- 
steady, zero-viscosity, irrotational theory of airfoils, which was already available at 
the beginning of the twentieth century. 

Classical airfoil theory explained satisfactorily the balance of forces in a wing in 
steady flight. In short, the force that air exerts on the wing is divided into two 
standard components: the lift (vertical force) and the drag (horizontal force), where 
horizontal means the direction of steady motion. In steady flight, these forces are 
balanced by the weight and the propulsion force. The theory predicts that the lift 
and the drag are proportional to the circulation of air velocity around the airfoil, 
and it was in agreement with experiments, see [1] for details. Trouble occurs when 
one wants to change the lift, as one should do when attempting to take off or land 
in a fixed wing aircraft. A theorem, due to Lord Kelvin, states that the circulation 
around a material curve, such as the boundary of an airfoil, is a constant of motion 
in ideal (i.e. non-viscous) flow - or maybe, nearly constant in slightly viscous 
flow. So, changing propeller speed and moving control surfaces does not change 
the circulation. Since airplanes start out at rest, with zero circulation around the 
wings, no airplane could, on theoretical ground, develop a lift, and therefore fly. 
Something was clearly wrong with the theory. 

The correction was due to the young theoretical mechanician Ludwig Prandtl 
(1875 - 1953), who published a short paper in the Proceedings of the Third ICM 
(Heidelberg 1904) whose German title roughly translates as "fluid flow in very 
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little friction". In this article, Prandtl established a perfectly satisfactory, and 
revolutionary, explanation of the following observation: 

(O) The interaction of incompressible flow with a material boundary is com- 
pletely different if the flow has very small viscosity or none at all. 

This observation, the associated explanation, called boundary layer theory and 
some of what mathematicians made of this subject in the following century and a 
bit, make up the subject of these lectures. For a thorough account of the develop- 
ment and understanding of the physics of boundary layers, we would like to refer 
the reader to the classical book 29J. It can be argued that this short paper by 
Prandtl marks the birth of modern applied mathematics. 

The theory of boundary layers is a cornerstone of modern fluid mechanics, but, 
as in much of this field, it lacks a rational framework, i.e. a rigorously established 
connection with first principles. Although substantial mathematical work has been 
done in this direction, some basic questions remains unanswered. The purpose of 
these lectures is to probe the boundaries in the mathematical understanding of 
the interaction between nearly ideal flow and solid objects, perhaps to bring what 
is not known about this question more sharply into focus. The choice of material 
covered is strongly slanted towards recent work by the author and his collaborators, 
and it includes detailed consideration of Kato's criterion for the vanishing viscosity 
limit in a bounded domain, a long discussion on the vanishing viscosity limit for 
incompressible flows with Navier boundary condition and the detailed behavior of 
circularly symmetric flow inside a rotating cylinder. The choice of working with 
two dimensional flow is both a reasonable pedagogical choice and a comfort zone 
for the author - in the issue of boundary layers, the sharp distinction in behavior 
between 2D and 3D flows is not yet apparent, and much of the work we will discuss 
here generalizes readily to 3D. Finally, we mention that these notes are written 
thinking of a reasonably mature audience - we assume, not only familiarity with 
standard PDE theory, but some familiarity with the basics of mathematical fluid 
dynamics as well. 

The remainder of these notes is divided in seven sections as follows: Section 2 
contains a derivation of Prandtl's equation; Section 3 contains a broad overview 
of rigorous results on Prandtl's equation, including some of O. Olcinik's work, and 
more recent progress; Section 4 introduces and proves Kato's criterion and some 
related results; Section 5 is concerned with vanishing viscosity under Navier friction 
conditions and a proof of L p vorticity estimates in this case; Section 6 contains an 
exposition on a rigorous method to treat boundary layer expansions based on ideas 
of geometric optics, applied to the Navier boundary condition; Section 7 explores 
a nearly explicit example of the behavior of the boundary layer for the no-slip 
condition; Section 8 contains some conclusions and open problems. 



In this Section, we present an asymptotic derivation of Prandtl's boundary layer 
equation. Our point of departure is the Navier-Stokes equations, which are an 
expression of Newton's second law applied to the motion of a fluid, subject to an 
incompressibility constraint. We write 



2. Prandtl's theory 



(1) 



dtu + u ■ Vu 
div u = 0, 
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where u is the fluid velocity, p is the scalar pressure and /i is the kinematic viscosity 
of the fluid. We have assumed that mass units have been chosen so that fluid density 
is one. 

For this derivation, and for most of the discussion in these lectures, we will 
assume that we are discussing a two-dimensional fluid occupying a half plane H = 
{x2 > 0}. Two dimensional fluids really occur either in a computer simulation, or as 
three dimensional fluids which are translation-invariant along some direction (often 
an unstable state of affairs) or as an approximate model for thin fluid layers. Much 
of the discussion extends very naturally to a fairly general fluid domain in three 
dimensions, but we will stay with the simplest possible situation for pedagogical 
reasons. 

As we all know, problem {T]) requires a boundary condition at {x 2 — 0}, and the 
condition usually deemed appropriate is the no slip condition u(xi,0,t) = 0. This 
condition expresses an assumption that viscous fluids adhere to material objects, 
something that is neither physically nor mathematically obvious and was subject 
for heated debate until the mid nineteenth century, when it became clear that it 
gave good agreement with experiments. 

Ideal, or inviscid flow is represented by solutions of Euler's equations, which 
is system ([T]) with /i = 0. For ideal flow, the correct boundary condition is the 
non-penetration condition U2{x\,0,t) = 0. The Navier-Stokes system is a singular 
perturbation of the Euler system, because the small constant fi appears in front of 
the highest order term of ([I]) . One consequence of this singular perturbation is the 
disparity in boundary conditions between Euler and Navier-Stokes flows - namely 
that Mi at the boundary goes from being identically zero for any positive viscosity 
to some (in principle) nonzero function when /i = 0. This disparity is the root 
cause of the boundary layer trouble. 

Our objective here is to derive Prandtl's boundary layer equations. This is a 
way to quantify nearly ideal fluid behavior near the boundary by means of an 
appropriate set of limit equations. Let us begin with a simplifying assumption: we 
assume that the disparity between ideal and viscous flow is concentrated on a thin 
layer near x 2 = 0. 

Our next step is to non-dimcnsionalize equation ([I]) , using the time scale T, the 
length scale L for horizontal lengths, and a reference vertical length scale h. We 
introduce the non-dimensional constant v — fiT/L , which is a measurement of the 
quotient of viscous by inertial forces in our flow, and measures in a physically ap- 
propriate manner how far from ideal our flow really is. The non-dimensionalization 
procedure simply means introducing new variables 

u 2 {y,s) 

and 

p(y, s) 

which results in the system 



T 

= ju {Lyi,hy 2 ,Ts), 
T 

T 2 
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(2) 



dgU 1 + u ■ VyM 1 = —d Vl p + v 



I? 



dl u 1 + ^dlu 1 



L' 2 



d s u 2 + u ■ V y u 2 = --^d y . 2 p + V 



d 2 u 2 + —d 2 u 2 

Vl h V2 



divj,w 



0. 



We introduce e = h/L, which is assumed to be a small, non-dimensional, param- 
eter because we are focusing in a thin layer. We also want to consider v small. A 
key issue in boundary layer theory is that the magnitude of the small parameters e 
and v are naturally related. Indeed, if v <C e 2 then matched asymptotics indicates 
that, to leading order, u satisfies the Euler system (with pressure independent of 
the vertical variable) and with no-slip conditions at the boundary. These boundary 
conditions are inconsistent with the fact that the Euler system is of first order. On 
the other hand, if v 3> e then, to leading order, u is such that d 2 2 u = 0, which, 
together with the no-slip boundary condition implies that u = c(yi)y2, for some 
vector c. Now, if u is to represent the behavior of the flow in a thin layer near 
the boundary, then the velocity u should match the inviscid velocity in the limit 
1/2 — > oo, and not just blow-up. The only regime that appears to yield a consistent 
asymptotic regime is 



(3) 



v/e 2 =0(1). 



From another perspective, condition ([3]) highlights the region near the boundary 
where the vertical viscous stress balances the inertial terms in the Navier-Stokes 
system. Assuming v = e 2 and implementing matching asymptotics for v small 
we obtain the following system for the leading order approximation, denoted v — 

{v\v 2 ), 



(4) 



dsV 1 - 

d y2 q = 

db/yV = 0. 



These are the unsteady Prandtl equations for the boundary layer profile v. They 
represent the behavior of the flow near the boundary. To obtain a complete problem, 
these equations must be supplemented with boundary conditions. First we impose 
the no-slip condition 

v = at ?/2 = 0. 

An additional condition must be imposed in order to capture the assumption that 
far from the boundary layer the small viscosity Navier-Stokes solutions match the 
Euler solutions. Let u v a family of solutions of the non-dimensional Navier-Stokes 
equations with non-dimensional viscosity v and u E be a solution of the incompress- 
ible Euler equations. For example, can assume that both the family u v and u E are 
defined by solving the Navier-Stokes equations and the Euler equations with the 
same initial data u u {x,Q) = u E (x,0) = uq(x). 

Going back to the Prandtl system, we expect that v(yi,y2,t) — > u E (yi,0,t) = 
U(yi,t), as ?/2 — ► oo. Let p E be the pressure associated with the Euler solution u E . 
Since q does not depend on j/2, looking at j/2 — > oo makes it also natural to assume 
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that q(yi, t) — p E (yi, 0, t). If we look at the Euler equations and evaluate them at 
x 2 = we obtain the following relation: 

U t + UU yi = -pE = -q yi , 

which is called Bernoulli's Law. This means that, for the Prandtl equation the 
condition at infinity U determines, up to an irrelevant constant, the pressure q. 
With this construction, we hope that, when v is small, 

tx\ ..u( r _ * _ J (v 1 ,ev 2 )(xi,x 2 /e,t) for x 2 < \(v) , 

(5) u [x u x 2 , t) - | u E {xuX2ft) for X2 > x{y) + oil), 

where A(^) is any cutoff distance such that e <C \{v) <C 1. 

In addition to the time-dependent Prandtl equation, this derivation also yield 
the steady Prandtl equation, given by 

!v ■ Vyv 1 = -d yi q + d^v 1 
d y2 q = 
dWyV = 0. 

The typical problem associated with this equation is a quarter-plane BVP, where a 
profile v(0, y 2 ) is given and one attempts to find v(yi,y 2 ) for y\ > as the induced 
boundary layer profile over a half-plane plate. 

The derivation above is a nice example of multiscale asymptotic analysis, and 
from the complicated issue surrounding the interaction of nearly inviscid flow with 
material boundaries it derives a new equation, ([5]), and a simplified asymptotic 
model for the behavior of Navier-Stokes solutions near a material boundary, given 
by ©. The key issue that such a model raises is its validity (mathematical) and 
applicability (physical). 

This model has been found useful in applications, specially where it concerns 
laminar boundary layers, and when its usefulness begins to break down, suitable ex- 
tensions of the model have been obtained, notably the so called "triple deck" expan- 
sions, where an intermediate thin layer is added between the viscosity-dominated 
internal layer and the free irrotational Euler flow. In this intermediate layer, the 
flow is ideal, but not necessarily irrotational. One important situation where the 
boundary layer ansatz breaks down is when boundary layer separation occurs. Re- 
call that one of the assumptions in deriving the Prandtl equation was that the 
disparity between ideal and viscous flow be concentrates in a thin layer near the 
boundary. It is quite common, even well within laminar flow regimes, that the 
boundary layer detaches itself from the boundary and affects the bulk of the flow. 
In that case, Prandtl's theory and its extensions break down as models. In the 
next section we will consider what is known regarding the rigorous validation of the 
asymptotic approximation described here. 

3. PRANDTL'S EQUATION 

The purpose of this Section is to present the known theory for Prandtl's equation, 
without much detail. The first question one must address regarding an approximate 
model is: can I solve it? The initial and boundary value problems for the Euler 
and Navier-Stokes equations are well-posed, globally in the case of the half-plane 
with reasonable initial data. So the issue of whether Euler + Prandtl is a good 
approximation for Navier-Stokes with v small, in the sense discussed in the previous 
section, depends first on understanding the well-posedness for Prandtl's equation. 
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The mathematical theory of the Prandtl equation only got started in the sixties, 
by O. Oleinik. Over the years, Oleinik and her group made a large number of 
contributions to the theory of Prandtl's equation and its many variants, collected 
and explained in the book |27j . For the present discussion, we would like to focus 
on one specific result, that first appeared in [26]. We also refer the reader to the 
survey [5], for the discussion of Oleinik's result and its relation with the blow-up 
result of W. E and B. Engquist, [6]. 

Let v — (v 1 , v 2 ) be a solution of the IB VP for Prandtl's equation, which we write 

as 



(7) 



dtv 1 + v ■ Vv 1 = -d Xl q + d 2 2 v Y 
divv = 

v{x\ , 0, t) = and lim^^oo v 1 {x, t) — U {x\ , i) 
v(x, 0) = v (x), 

Where —d Xl q — Ut + UU Xl . The result we wish to discuss is the following 

Theorem 1. (Oleinik 1967) Assume that both U and Vq are positive and that, in 
addition, 8 X2 Vq > 0. Then there exists a unique global strong solution of ([7]). 

We will not present a proof of this result, but we will discuss a key part of the 
proof, which is the recasting of this problem as a scalar, degenerate parabolic scalar 
equation, using Crocco's transformation. 

We begin by taking the derivative of Prandtl's equation with respect to xi and 
introduce the new dependent variable u(x,t) = d X2 v . System ([7]) is equivalent to 
the following IB VP: 



(8) 



d t u) + v ■ X7uj = dl 2 uj 
v = K[oj] 

(d X2 Lj)(x 1 ,0,t) = d Xl q and lim l2 ^ 00 w(i,t) = 
oj(x,0) = d X2 v (x), 

where the vector operator K reconstructs (v , v 2 ) by first integrating in the vertical 
variable to obtain v\ and then using the divergence-free condition and integrating 
again in the vertical variable to obtain V2- The new equation ||HJ| is, in a sense, the 
vorticity formulation of problem (|7|. 

We assume that the solution v — v(x, t) we seek satisfies the condition d X2 v 1 (x, t) > 
for all x £ H, t > 0. In particular, this means that, for each fixed (x\, t), the map 
x-2 I— > v (xi, X2, t) is invertible. Let us denote this inverse by h — h(xi,£,t). In 
other words, we have 

(9) v 1 {x x ,h(xi,£„t),t) =£, for all £ > 0. 

The Crocco's transform consists of introducing the new dependent variable: 

(10) W = W{xx,£, t) = lj( Xi , h(x u ^ t),t), 
We verify that W is a solution of the following IB VP: 

r 8 t W + iW Xl - (d Xl q)W i = W 2 d 2 W 

(11) <^ WW £ = q Xl for e = 

[ W{x,0) = (d X2 v Q )(x 1 ,h(x ll ^0)) 

Indeed, we can compute directly to obtain: 
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d t W = — t-u^+ojt; d Xl W 




,2 
,3 



Substituting the corresponding equalities above into (fTTj) . and using (fTUJ) , © 
and the evolution equations in ([7]) and © verifies the evolution equation in (fTTj) . 
In addition, we can check directly that 

u> X2 (xi,x 2 ,t) = W(xi,g~ 1 (xi,X2,t),t)W^(xi,g~ 1 (xi,x 2 ,t),t), 

which, together with the boundary condition in (jSJ) gives the boundary condition 
in (fTT|) . Problem pip is a scalar, degenerate parabolic equation, which is amenable 
to fairly standard treatment, using fixed point methods, and satisfies a comparison 
principle. In particular, the sign of W is retained in the evolution, and therefore, 
the monotonicity condition on v , necessary for the validity of Crocco's transform, 
is retained as well. 

Finally, once a solution W is obtained for problem (JTTJ , one must reconstruct a 
solution to the original problem. Recall that 



which allows the reconstruction of u from W by means of (fTTjj) . The interested 
reader may prove, as an exercise, that such an u) is, in fact, a solution of (|5]l. 

This result, and others proved by Oleinik and her group, give useful, rigorously 
established descriptions of the vanishing viscosity asymptotics, but depend, to a 
greater or lesser extent, on monotonicity conditions such as 8 X2 Vq > 0. As we have 
seen, the monotonicity assumption is needed for the validity of the Crocco's trans- 
formation, but this assumption might just be a feature of the method, rather than 
an essential limitation of the theory. In 1997, E and Engquist produced a coun- 
terexample which showed that Prandtl's equation develops finite-time singularities 
if the monotonicity condition is not imposed, see [6]. In fact, E and Engquist's ex- 
ample suggests that the role of the monotonicity assumption is to prevent boundary 
layer separation, a phenomenon that actually occurs in real flows and corresponds 
to a breakdown of the Prandtl ansatz. 

An alternative to the half-plane analysis described above is to study well-posedness 
of Prandtl's equation in bounded intervals in x%, where the horizontal velocity of 
the boundary layer profile is specified in one side of the interval, and the length of 
the interval or the time horizon of the analysis are chosen small enough to prevent 
separation. Such a result was first proved by Oleinik in |25j . Recently, Z. Xin and 
L. Zhang improved Oleinik's result, showing existence of a global (in time) weak 
solution for Prandtl's equation on a finite horizontal interval, if the pressure is fa- 
vorable, i.e., q Xl < 0, see |34j and |35j . This condition is also known to discourage 
boundary layer separation. 



v 1 {x,g(x 1) £,t) ) t) = £, 
and therefore, differentiating this identity with respect to £ gives 

dg 

u>{xx,g{xi,£,t),t) — {xi,£,t) = 1. 
Recalling pTj)) . this implies that 
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A different approach to the theory of Prandtl's equation was taken, initially by 
A. Asano, in a couple of unpublished manuscripts, and later by R. Caflisch and M. 
Sammartino, in a pair of articles, see 28J, recently further improved by Lombardo, 
Cannone and Sammartino in [17] . The basic idea is that, without the monotonicity 
condition, or something analogous to it, one expects the initial-boundary value 
problem ([7|) to be ill-posed. As a result, it becomes natural to look for local (in 
time) solutions for Prandtl's equation in analytic function spaces, using results 
of Cauchy-Kowalewska type. The main results in [28) were well-posedness of the 
problem ([7]) if the data vq and U are analytic, and compatible. In [17] the analyticity 
requirement on vq was imposed only in the horizontal variables. 

Of course, the well-posedness in analytic spaces, and the blow-up example by E 
and Engquist does not prove that ([7]) is ill-posed, which at this time remains an 
interesting open problem. 

To conclude this section, it would make sense to mention the contribution of 
E. Grenier, which he describes as nonlinear instability of the Prandtl boundary 
layer. His result is not about Prandtl's equation per se, but about the vanishing 
viscosity limit of the Navier-Stokes equations. His result can be interpreted as 
mathematical evidence that the Prandtl ansatz is not always valid for solutions of 
the Navier-Stokes system in the half-plane with small viscosity, see [7J. In other 
words, although the theory of Prandtl's equation is relevant for understanding the 
vanishing viscosity limit, there is more to the original observation (O) than Prandtl's 
original explanation for it. 

4. Kato's condition 

In this section we move away from Prandtl's equation, and we begin a study of 
the vanishing viscosity limit from a broader point of view. Our first observation 
should be that, even in the absence of boundaries, all the mysteries of turbulence 
lurk in the background of the vanishing viscosity limit, see for example [18 for a 
small part of this story. However, under moderate regularity assumptions, for exam- 
ple, if the initial vorticity is bounded, explicit estimates for the difference between 
Euler and small viscosity Navier-Stokes solutions are known, see [2]. Also, and this 
distinction is a key point here, for very irregular flow, we still have the existence 
of subsequences of solutions of small viscosity Navier-Stokes converging to weak 
solutions of the Euler equations, up to and including initial vorticities which are 
measures, see [23l [22] , but then no estimates on the difference are provided, or ex- 
pected. Basically, in the absence of boundaries, as long as the underlying ideal flow 
has enough regularity so that uniqueness of weak solutions to the Euler equations 
is known, we have actual convergence of the vanishing viscosity limit. Furthermore, 
as long as existence of weak solutions is known we also have compactness of the 
vanishing viscosity sequence and weak continuity of the Euler/Navier-Stokes non- 
linearity. Nonuniqueness of weak solutions for Euler equations is also known, see 
the remarkable paper [4], and references therein, for the current knowledge on this 
nonuniqueness, but the behavior of the vanishing viscosity limit for these examples 
is a very interesting open problem. 

As soon as we consider flows in the presence of boundaries, the story changes 
quite dramatically. Very little is actually known mathematically, and this very little 
is precisely the object under discussion in these notes. Physically, boundaries are 
the most natural source of small scales in incompressible flows, precisely through the 
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boundary layer mechanism, and these small scales are the source of the irregularities 
that justify considering irregular 2D flows in the first place. The point of departure 
in our discussion will be a classical open problem, which we formulate below. 

d Layer Problem: Let u v be a sequence of solutions of the incompressible 
Navier-Stokes equations in two space dimensions, in a smooth bounded domain Q, 
satisfying the no-slip boundary condition on dQ, with initial data Uq, bounded in 
L 2 . Is there a subsequence u Uk converging weakly in L 2 to a vector field u, which is 
a weak solution of the incompressible Euler equations in f2 with some initial data 
Mo = limuQ fe ? 

This problem is open even if luq = ujq G C£°(fi), with luq = curl uq. Let 
us focus, for simplicity, in this case. Clearly, the Navier-Stokes equations have a 
unique smooth solution u v with initial data uq, and the Euler equations also have 
a unique smooth solution u with the same initial data. We will see that there are 
examples where u v — > u in L 2 , but the answer to the problem above may be positive 
even when u v does not converge to u, because there may be weak solutions of the 
incompressible Euler equations with initial velocity uq which are not u. 

In 1984, T. Kato wrote a short note where he proved a sharp criterion for the 
convergence of u v to u, see [11]. The observation by Kato is remarkable for at least 
two reasons. First, as we shall see, it is very natural from the analytical point of 
view. Second, it places the condition for convergence on the behavior of the small 
viscosity sequence at a distance 0(y) of the boundary of 17, hence in a much smaller 
region than what is the natural domain of the boundary layer. Next, we state and 
prove a simple version of Kato's criterion. 

We focus in the case fl = {\x\ < 1} in K 2 . Let lu e C£°(fl) and u a = K[uj ], 
where K is the Biot-Savart operator in the unit disk. Let u v be the unique classical 
solution of the Navier-Stokes equation in with no-slip boundary condition and 
initial velocity uq, and u be the unique smooth solution of the Euler equations with 
u ■ x = for |x| = 1 and initial velocity Uq. 

Theorem 2. (Kato 1984) Fix T > 0. There exists a constant c > such that 
u v -> u strongly m L°°((Q, T); L 2 {VL)) if and only if v \\Vu v \-,t)\\\ 2{T ,dt -> 
as v — > 0, where T cu = {1 — cv < \x\ < 1}. 

Proof. First consider the energy identities for both u v and u. We have, for each 
t > 0, 



ll u 0»*)lll=(fi) = IKHl2 ( fi). 
Therefore, if u u -> u strongly in L°° ( (0, T) ; L 2 (O) ) , then 1 1 u v ( • , t ) | \\ 2 (n) -> 1 1 u ( • , t) \ 




and 



for almost all time, and therefore 
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for each t > 0, not almost everywhere anymore since the integral in time is increas- 
ing in time, and therefore, 

v I I \Vu"\ 2 dxdt -» 0, 
Jo ir„ 

as we wished. 

To prove the other implication, fix e > and let <\f € C£°(Q) be such that 
(f> £ {x) = ip E (\x\), with ip e (x) = 1 for |ar| < 1 - e, ^ e (x) = for 1 - e/2 < |x| < f , 
and (f e decreases monotonically from f to 0. Let u — curl u be the vorticity and 
V> be the stream function associated with the Euler flow u. Define 

u e = V^V) = (-5 X2 (^v),s xl (^V))- 

Let u £ = u — m £ = V ± ((f — (f> £ )ip). The stream function ^ vanishes at |x| = 1, 
and it can be assumed to be uniformly bounded in C fc , for any k = 1, 2, . . ., so that 
we can easily obtain the following estimates on v £ : 

(12) \\v e \\L°°((a,T);L 2 (n)) < Ce 1 / 2 

(13) ||St^|Ui( ( o,T);L2 ( a)) ^Ce 1 / 2 

(14) ||Vw e || L oc (( o,D;L»(n)) < ^" 1/2 

(15) |KIU°°((0,T)xfi) < C 

(16) l|Vw £ || L oo ((0 iT)xQ) < Ce" 1 

In addition, we require certain estimates on u", uniform in i/, which we collect 
below 

(17) |K||loo( ( o,T);L2(Q)) < C 

(18) "l|V«"||ioc (( o,r);^(n)) < C. 
By using Cauchy-Schwarz in time, we also have 

(19) ^ /2 ||V^|U 1(( 0,T);^(O)) < CT 1 ' 2 [v ||V^|I^(Q)^ < C 

Finally, a version of Poincare's Inequality, which reads 

(20) KIU'(r.) < CeHV^IU^r.) 

Now we estimate, omitting the explicit dependence of v on e: 

IK - u\\l Hn) = \\u»\\l Hn) + |M|| 2(n) - 2(u\u) 
< 2\\u \\ 2 L2{n) -2(u»,u-v)-2(u»,v). 
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We have that 

\{u V ,v)\ < ||u'lL~((0,T):L2(n))|M| L oo(( 0iT ). L 2( O )) < Ce 1/2 . 

And therefore, taking e = Cu, 

(21) \\u» - u\\ 2 LHn) < 2\\u \\ 2 LHn) - 2(u", u-v)+ 1(1). 

We multiply the Navier-Stokes equation by u — v , integrate in space and time, 
and integrate by parts to obtain 



(u v , u — v) — (uq, uq — v) = / (u u , u v ■ V(u — v)) — j/(Vti", V(u — v )) + (u v , dt (u— v))dt, 

Jo 

which, together with (|2~lj) implies 
(22) 

ll«"-»lli»(n) < / [-2(u" ^ '-V{u-v))+2v(Vu v ,V(u-v))-2(u v \d t (u-v))]dt+o(l). 



We write 



and we have 



u v , dt(u — v))dt = 
o Jo 



(u l/ ,u-Wu)dt + / (u u ,d t v)dt, 



(m", d t v)dt 



< / || u || L 2||a t t,|| L2 <cv 1/2 , 

Jo 



hence. 



(23) 

IK-^lli=(n) < / [-2(u" ,u v ■V{u-v))+2v(Vu" , V(u-v))+2{u v ,u-Vu)]dt+o{l). 
Jo 

Note that 

{(u u -u),(u u -u) ■ Vit) = (u v ,u v ■ Vu) - (u v ,u- Vu), 
and therefore, ((23)) becomes 

(24) 

\\u v -u\\ L 2 {n) <2 [ [~{(u v -u),{u v -u)-\7u)+v{\7u 1 ' ,V{u-v))+{u v ' ,u v -Vv)}dt+o{l). 
Jo 



We analyze each term in the right hand side of (|24|) separately. First 



{(u v - u), [u v -u) -Vu)dt 



<l \\u v -uf mm \Wu\dt<C I \\u v -uf mm dt. 



Second, 



v \ {Vu v ,V{u-v))dt 



< 



v \ (Vu",Vu)dt 



v \ (Vu v ,Vv)dt 



< 



v / l|Vu ly || L 2(fi)l!V'u|| L 2 (n) cte + v I ||\7u ,y || i 2 (rciy) ||\7u|| L 2 (rci/) 
Jo Jo 

< C^||Vu i/ || l i(( 0i t);L2(o)) + C , ^ 1/2 ||Vu l/ || L i(( 0: T);L2(r c „)) 
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< Cv^H^Vu^U^Ty^^)) + o(l) = Ct 1 ' 2 Lj \\Vu"\\l, [T ^ 



1/2 

= 0(1), 



by Kato's criterion. 
Third, 



t 

(w , w 



t 

(v, vr 



o 





t rt 



< 2 ||w|U«>((o,T)xn) / l|w l/ ||L2(r c „)l|Vu' y || L 2 (rcw) ^ < Cv / HViil^r^cft = o(l), 

Jo Jo 
where we used Poincare's Inequality and Kato's criterion in the last line. 

Therefore, we get back to (j24)l and use the estimates obtained above to conclude 
that 



ft 



% m <C I ||«"-u||i, (n) £ft + 0(l), 

'o 



which, by Gronwall, concludes the proof. 



□ 



A number of observations are in order. First, the proof above can be adapted to 
general bounded domains in R™, to Leray solutions of the Navier-Stokes equations 
and to problems with forcing and different initial data, as long as the forcing and the 
initial data converge to the corresponding ones when viscosity vanishes. The two 
equivalent conditions in the statement of the Theorem are also equivalent to weak 
convergence of u v to u in L 2 , pointwise in time. Clearly the strong convergence 
implies the weak convergence, and we weak convergence implies that 

HVuli^njdt - 0, asz/^O, 

o 

which can be seen using the energy identities and the weak lower semicontinuity of 
the L 2 norm. All these observations were contained in the original work by Kato, 

HQ- 

There have been several results modifying, extending and improving Kato's cri- 
terion in the literature: work by R. Temam and Xiaoming Wang, |30[ 31 , where 
Kato's condition is replaced by a condition on integrability of pressure near the 
boundary, work by Xiaoming Wang, [33] where the criterion is cast in terms of 
integrability of derivatives of tangential components of velocity, and work by James 
Kelliher, [13], where the integrability condition is placed on vorticity. 

As we will see in the last lecture, there are some (trivial) examples where Kato's 
criterion holds true, but this is not known in general, and it would not be surprising 
if it turns out to be false most of the time. Note that, if Kato's criterion fails, this 
implies non-vanishing energy dissipation occurring in a thin neighborhood around 
the boundary, as viscosity vanishes. 

Given the disparity of the understanding of the vanishing viscosity limit in the 
presence or absence of boundaries, it becomes natural to ask how much boundary it 
is required to obstruct the proof of convergence of the vanishing viscosity limit. The 
author has been involved with two recent results in this direction. First, with D. 
Iftimie and H. Nussenzveig Lopes, we have proved that if u^' 5 is the solution of the 
Navier-Stokes equations in the exterior of the domain efl, with no-slip boundary 
conditions and fixed initial vorticity and u is the solution of the Euler equations in 
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the full space, with the same initial vorticity then there exists a constant c > such 
that, if e < cis, then u u ' £ — > u strongly in L 2 , see [8]. Second, with J. Kcllihcr and H. 
Nussenzveig Lopes, we have proved that if u u ' R is the solution of the Navier-Stokes 
equations on a domain which contains the ball of radius R, with fixed initial 
vorticity, and again u is the solution of the Euler equation with the same initial 
vorticity. We have proved that, as long R — > oo as v — > 0, then u v ' R converges 
strongly to u as v — > 0, see [14]. Both proofs follow the basic idea and the structure 
of the proof of Kato's criterion as we presented it here. 

5. Vanishing viscosity with Navier friction d condition 

The main purpose of this Section is to highlight the role of vorticity production 
by the interaction of the flow with the boundary in the difficulty surrounding the 
Boundary Layer Problem. In particular, we briefly discuss the vanishing viscosity 
limit for flows in the full space, and also for flows in a bounded domain satisfying 
Navier friction condition at the boundary. Let us begin with the full space case. 

A basic ingredient for the proof of compactness of the vanishing viscosity limit 
in the style of [23] is the set of estimates obtained through the vorticity formulation 
of the Navier-Stokes equations. Vorticity play a fundamental role in the analysis 
of incompressible fluid flow. This is specially true in two dimensions, because 
vorticity is conserved along particle trajectories for ideal flow, and it satisfies a nice 
convection-diffusion equation for viscous flow, namely: 



where K is the Biot-Savart law associated with the flow domain. In the absence 
of material boundaries, for example, in the case of the full plane, equation (|25j) 
implies a priori L p estimates for vorticity for any 1 < p < oo, as long as the initial 
vorticity is in L p . The Biot-Savart operator is a pseudo-differential operator of 
order —1, which implies compactness of the sequence of velocities in L q , for any 
q < p* , p* — 2p/(2 — p) (assuming 1 < p < 2), and the nonlinearity of the Euler 
equations being quadratic, all that is required is compactness in L 2 in order to 
prove that the flows obtained through the limit of vanishing viscosity satisfy Euler 
equation. The case p = 1, which is what is involved in the vortex sheet initial data 
problem, is critical for the argument above, and requires a more involved analysis. 

To highlight the importance of the added compactness provided by vorticity 
estimates, we consider briefly the situation without it. In that case, the only a 
priori estimate available comes from the energy estimate, which gives solely an 
L 2 -bound on velocity, independent of viscosity. By Banach-Alaoglu, an uniform 
L 2 -bound can provide us with a sequence of approximations which converge weakly 
in I? to a "candidate" for limit solution u. To prove that u really is a solution of 
Euler 's equation, one must pass to the limit in the nonlinearity M n -Vn", which, after 
integration by parts on a weak formulation, can be written in terms of quadratic 
expressions in the components of velocity. Now, the point is: quadratic expressions 
are not continuous with respect to the weak topology of L 2 . Two classical examples 
(in one space dimension) are: 



(25) 





(concentration), 
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and 



g n = sin(nx) (oscillation). 



Both these sequences converge weakly to zero in L 2 , but their squares converge 
to something else - we leave the details as an exercise to the reader. The bottom 
line is this: strong convergence in L 2 is required to pass to the limit in quadratic 
expressions, and this strong convergence has to come from some additional a priori 
estimate. 

The presence of vorticity estimates breaks down in the presence of material 
boundaries because, although vorticity still satisfies equation (|25p in a domain with 
boundary, there is no natural boundary condition for vorticity to complete equation 
(|25|) . beyond the nonlocal boundary condition K [us] ■ r = (only the vanishing of 
the tangential component of velocity is required to vanish because the vanishing of 
the normal component of velocity is implicit in the Biot-Savart law). In particular, 
one may argue that it is the lack of control on the production of vorticity at the 
boundary which obstructs the solution of the d Layer Problem. 

Mathematicians are not discouraged by difficulties - if we cannot solve a problem, 
we find another problem nearby which we can actually solve. In the book 16J, J.-L. 
Lions used as an illustration the following problem: 



which is usually called the "free boundary problem" for the Navier Stokes equa- 
tions in the planar domain f2. The standard theory for the Navier-Stokes equations 
with no-slip boundary condition can be easily adapted to this problem, but even 
more, the vanishing viscosity limit is quite well behaved in this case. Indeed, the 
vorticity equation ([23)) can be closed by adding the Dirichlet boundary condition 
u> = at dfl and initial data. By multiplying the equation by a/'" 1 , integrating, 
and performing the usual integration by parts, one quickly arrives at the conclusion 
that the L p -norms of vorticity are bounded independently of v, and the compactness 
argument outlined in the beginning of the section leads to an affirmative answer to 
this variant of the d Layer Problem. The homogeneous Dirichlet condition is dissi- 
pative for the heat equation; in particular, it means that the L p norms of vorticity 
must decay in time. 

There is a generalization of the free boundary condition called Navier friction 
condition which is both physically and mathematically interesting. The physically 
correct definition involves the rate-of-strain matrix of the flow at the boundary 
but, in two space dimensions, one can prove that this boundary condition can be 
written as u> = (2k — a)u ■ t, where k is the curvature of the boundary, a is the 
friction coefficient and r is the counterclockwise unit tangent vector to d£l, if f2 is 
a bounded domain with connected smooth boundary. We refer the reader to [3] for 
a full discussion. This condition, also called Navier "slip" condition, still allows for 
the creation of vorticity at the boundary, but in a way that can be controlled. 

In [3], T. Clopeau, A. Mikelic and R. Robert studied the initial-boundary value 
problem for the Navier-Stokes equations in 2D, with Navier friction condition. They 
proved global well-posedness for this problem for fixed viscosity, and they proved 
the convergence of the vanishing viscosity limit to the (unique) solution of the Euler 



(26) 
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equations if the initial vorticity is bounded. An estimate on the rate of convergence 
was obtained by J. Kellihcr in [12] . 

In [3T] the author, together with H. Nussenzveig Lopes and G. Planas, extended 
the vanishing viscosity result by considering flows with initial vorticity in L p , p > 2. 
The key point of this result is the uniform L p estimate on vorticity, which is what 
we will do next. Let be a bounded, connected and simply-connected smooth 
domain in M. 2 and a S C°°(cMl) . We consider to" = oj v (x,t) to be the unique 
solution of the IB VP: 



(27) 



u v = K n [uu»] 

lo u (x, t) = (2k - a)(u(x, t) ■ t) for x G dQ 
to(x, 0) = loq 



We assume that ojq G H 1 n L°° is compatible, i.e. that ojq = (2k — a)(ifo[wo] ■ r) 
in dfl. Existence of a solution to problem (g7]| with oo v E C([0, T); H l (Cl)) (and 
u v G C([0, T); H 2 (Q) j) was proved in [3], and density of compatible vorticities in 
V with respect to the strong norm was proved in |21j . Our objective is to prove 
the following 

Theorem 3. (L., Nussenzveig Lopes and Planas, 2005) For each p > 2, there 
exists a constant C > 0, independent of v , such that 

\\u"(-,t)\\ L p {n) < C(||w |Up(fi) + \\u \\ L 2 (n )), 

with u = Kq, [lu ] . 

Proof. In this proof we will use some results of the theory of parabolic equations, 
which can be found in the book 15J. 

The proof involves applying a maximum principle to two auxiliary problems. 
First observe that u ■ r is H 2 , and therefore bounded. Set 

A = \\(2k - a)u ■ r|j L oo( a f2 X ( 0i T))- 
Consider the initial-boundary value problem for the Fokkcr-Planck equation: 



!u>t - vNJj + u ■ Vw = in tt x (0, T), 
2(-,0) = |wo| inO, 
w = A ondQx(0,T). 

This problem has a unique weak solution uj £ L 2 ((0,T); H 1 ^)). Then, lui = 
uj — to is a weak solution for the following initial-boundary value problem: 



{(wi)t - i/Awi + u ■ Vwi = inOx(0,T), 
wi(-,0) = wo-|wo| infi, 
u>i — (2k — a)u • t — A on 9fix (0T). 

The coefficients of the Fokker-Planck operator 9 t — vA + u ■ V are such that the 
maximum principle for weak solutions of parabolic equations is valid. Therefore, as 
u>i < on the parabolic boundary dil x (0, T) U fl x {t = 0}, it follows that wi < 
a.e. in f2 x [0, T). Analogously, we prove that u>2 = —to — lo is non-positive. We 
thus obtain 



16 



MILTON C. LOPES FILHO 



(30) M < ui a.e. in Q x [0,T). 

Moreover, since ujo is bounded, the maximum principle may also be applied to 
equation (gSJ) yielding that ui € L°°((0,T) x 0). 

Next we obtain an estimate for ui. Let cD = ui — A. This is a solution of the 
following problem: 

!Q t - vt±u + u ■ Vuj = infix(0,T), 
w(-,0) = |w |-A infi, 
£ = on<9Ox(0,T). 

We formally multiply (|31[) by cDp| p ~ 2 , where p > 2, we integrate by parts and 
use the incompressibility of the flow u to obtain: 

(32) ~ / \L>\P + (p-\)v [ |V£p| (p ~ 2)/2 2 dx = 0. 

In Jn 



Then, 



P dt , 

R-,*)||w(n) < |P(-,0)IU»(n) < IMIz»(n) +M^\ 1/p - 



Therefore, 

(33) PIIlp(o) < PIUp(o) + Airil 1 ^ < ||w ||^(n) + 2A|ft| 1 / p . 

This formal calculation can be made rigorous by using the weak formulation of 
(f3"Tj) given in [IS]. We begin by observing that cD t € £ 2 ((0, T); _ff _1 (f2)) and Q S 
L 2 ((0,T);i^(fi))nL°°((0,T)xQ). This implies that tD|w| p - 2 G L 2 ((0, T); iJ^(O)). 
Therefore we can multiply (f3"Tj) by a)|a3| p_2 if we understand the product with Q t 
and with AS as duality pairings. Finally, in order to justify (|32j) one still needs to 
approximate Q by suitable smooth functions and pass to the limit in each term of 
the weak formulation so as to obtain 

~ / \ujy = {Q t ,Q\Q\v- 2 ) 
pat J n 

= v(AQ,uj\Q\p- 2 } =-{p-l)v f ||V£|P| (p ~ 2)/ Tdx. 



Jn 

This can be easily accomplished using mollification in time together with the 
Dirichlet heat semigroup for f2, thus generating a family of smooth functions u) £ 
such that d t Q £ -> oJt strongly in L 2 ((0, T); H' 1 ^)), while uj £ \Q £ \ p - 2 Q\Q\ p - 2 
weakly in L 2 ((0, T); {ft)) and u e is uniformly bounded in fl x (0,T). 

Given (f33|) we now turn to the estimate of A. Using Sobolev imbedding and 
interpolating between W 1,p and L 2 , we find: 

K(-,i) ■ r|| Loo(an) < C\\u v {;t)\\c{h) < C\\u v (-MUa)\\ vV (-MwL in) 

<C\\u»(-MUn)\\" u (-ML;tn V 
where 9 = (p - 2)/(2p - 2). 

Let £ be an arbitrary positive number. We now use Young's inequality together 
with the fact that K and a are bounded to conclude that: 

(34) A < C , e ||« I/ ||£oo(( 0)T )._ E/ 2(n) + e||w !/ ||L~((o,T) ; Lp(o)) 
for some C £ > 0. Taking e small enough, from (|30 [) - (f34|) we obtain: 

(35) ||w 1, ||loo( 0jT . L p( )) < C(||wo||lp(o) + ||w' y |U=°(o,T;L 2 (n))) 
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for any p > 2, where C — C{p, Q, \\K>\\L°°(dCi)i ||c*||L°°(<9fi))- Finally, a standard energy 
estimate, yields ||u||l°°(o,T;L 2 (£-2)) < ||uo||i 2 (n)) thereby concluding the proof. 

□ 

Let us first observe that the " standard energy estimate" used as last step of the 
proof above is far from obvious. The difficulty lies in the integration by parts of the 
viscous term in the energy estimate. This can actually be done by using the correct 
formulation of the Navier friction condition, namely 2(DU) s n ■ t + era • r = 0, where 
(DU) S is the symmetric part of the Jacobian matrix DU. We leave the derivation 
of the energy estimate in this case as an exercise to the reader. 

The estimate in Theorem [3] can be extended to any vorticity in L p , p > 2 by 
density. It is not clear whether the restriction p > 2 is a consequence of the tech- 
nique used above, or it is an essential restriction in the estimate above. This is an 
interesting open problem. One physical justification for the Navier friction bound- 
ary condition is that is approximates the interaction of incompressible flow with 
a rough boundary. This contention was rigorously justified by Jager and Mikclic, 
who derived the Navier friction condition as the effective behavior associated with 
the homogenization of an oscillatory boundary, see |10j . 

There is no physical interpretation, at present, for the vanishing viscosity limit 
associated with Navier friction condition. It is, however, a very natural question 
from the mathematical point of view, since it interpolates naturally between the 
trivial case of the free boundary condition and the Boundary Layer Problem. The 
Navier friction condition still gives rise to a boundary layer, albeit in a less singular 
way than Prandtl's original theory. This is the subject of the next Section. 

6. Boundary layer for the Navier friction condition 

In this section we will explore a rigorous procedure to obtain a boundary layer 
equation, based on ideas from geometric optics, applied to the Navier friction con- 
dition. What we will present here is based on joint work of the author with D. 
Iftimie, H. Nussenzveig Lopes and F. Sueur. 

We begin with the i/-Navier-Stokes system 

' d t u" + u v ■ Vu" = -Vp" + Mil" 
div u» = 
(36) \ u%{t,xi,0) = 

u v (t,x 1) 0)=ju v 1 (t,x 1 ,0) 
u u (0,x) = u (x). 

We denote by u E , oj e , p E the solution of the initial-boundary value problem 
for the Euler equations in the half-plane, with initial data uq(x) and the standard 
no-penetration boundary condition. 

We set 



u = u + I q ) + vu R . 
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Additionally, we write 



u" = lj e + W{t, X!,X 2 /VP) + yfitj?, 
where W — W(t, Xi,y) = (d y v)(t, x±,y) and 4> v — \Jv curl u R . 

We are going to write an equation for v, derive an equation for W and then use 
the equations for W, lo v and lo e to deduce an equation for <j> v . Our final objective 
is to obtain an L p estimate for tp 1 " which is independent of v. 

We first write an evolution equation with Neumann boundary condition for the 
velocity profile of the boundary layer v: 

d t v + uf(t,X!,0)d Xl v + d Xl uf(t,X!,0)v + y{d X2 u E )(t 7 x 1 ,0)d y v = d 2 y v 
(37) I v(0,xi,y) = 

d y v(t,xi,0) = juf(t,xi,0) - cj s (i,xi,0). 

Now we differentiate equation (|3"T|) with respect to y and use the div-free condi- 
tion on u E to obtain an equation for the vorticity profile of the boundary layer W. 
We find: 



(38) 



' d t W + u E {t, Xl ,0)d Xl W + y{d X2 u E ){t, Xl , 0)d y W = d 2 y W 
W(0, Xl ,y)=0 



W(t,X!,0) ^juf(t,x 1 ,0)-uj E (t 7 x 1 ,0). 

The equations for v and W are linear, with smooth coefficients and independent 
of v. For the sake of these lectures, we assume well-posedness and exponential decay 
of the solutions v and W as y — > oo. To be more precise, we assume existence of a 
unique smooth solution v = v(t, x\,y) of (f57|) and consequently of W = W(t, x\, y), 
solution of such that: 

(1) W S C fc °° n LP, for some p > 2, 

(2) v e C b °° n LP and d xi v e LP, 

(3) d Xl W and d^W both in D>, 



(4) yd y W bounded and y 2 d y W e U>, 

(5) yd Xl W G LP. 

Next we write the equation for </)". First set z = X2 



We introduce 



and 



G" 



[uf{t, x u x 2 ) - u E (t, x u 0)](d Xl W){t, x u z) 



+ [u E (t,x 1 ,x 2 ) - x 2 (d X2 u E )(t,xi,0)} — (dyW)(t,X!,z), 



F v = V^Acj £ + y/v^ W) {t, x x , z) 



v(t, xi,z) 




Vw 1 



G" 



-(v(t,X!,z) + s/v(u v R ){)(d xl W){t,xx,z) - (u R ) 2 (d y W)(t,xi,z) -=. 

With this notation, the equation for <f> v becomes 
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' d t 4> v + u u ■ Vft - vAft = F v 

4> v (0,x) = 
t r(t,xi,0) = ^(v + ^(u R ) 1 )(t,x 1 ,0). 



(39) 



We break problem (|39j) into two problems, one with homogeneous boundary 
condition and nonzero forcing and another with nonhomogeneous boundary data 
and zero forcing. We write 4> u = <j)\ + 4>%, with 



(40) 



and 



' d t 4>\ + u v ■ Vft - vAft = 

<K(0,x)=0 
{ ft (t, Xl , 0) = 7 (v + v^K)i ) (*, xi , 0), 



' d t ft + u v ■ V<j% ~ v ^2 = FV 



(41) 



<%(o,x) = o 



The basic idea is to treat (|4"Tj) using the energy method and to treat (|4*D|) adapting 
the comparison principle argument used in 21j. In addition, it is necessary to obtain 
a priori estimates for the elliptic system which relates <j) v with u R . We write this 
system as: 

!div ^/vu R = d Xl v(t,Xi,xu/y/v) 
curl y/vu v R = ft 
K) 2 (t,xi,o) = o. 

What we end up with is the following estimates 

ll#IU»<c(i + ||#M, 

and 

^II^IU* < c{W\\ L , + 1) < C(i + n^iii, + ii^iu,), 

the proofs are the technical core of this work, but for time limitations we omit 
them. 

Just to conclude this story, the estimates above lead to 

j t Wi\\L, <C(1 + ||^IIl,). 

It follows from Gronwall's inequality that 11^2 II L p is bounded independently of v, 
for a fixed time interval [0, T]. Then, using the first estimate again concludes that 
there exists a constant C independent of v such that 



LP 



< c. 
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Remark: We observe that in [9] , Iftimie and Sueur proved that (ff is bounded in 
L 2 ((0, T)xfi), independently of v. We proved that there is a bound on L°°((0, T);LP), 
for 2 < p < oo. We do not expect to be able to extend this argument to p — oo. 
Note also that the analysis in [9] is valid for general smooth domains in two or 
three space dimensions, whereas the analysis outlined above only works in two 
dimensions. 

7. The rotating cylinder 

Our purpose in this section is to provide a fairly explicit illustration of the 
boundary layer phenomena. The main point of this illustration is to show, that in 
a favorable scenario, the small viscosity limit has vorticity behaving like a vortex 
sheet. This explains, in part, the difficulty involved in the Boundary Layer Problem, 
since vortex sheet regularity is critical for the weak continuity of the nonlinearity 
present in Euler and Navier-Stokes. This material is based on joint work of the 
author with A. Mazzucato, H. Nussenzveig Lopes and M. Taylor, and it is contained 
in two articles, pT9l [20] . 

The physical situation we wish to consider is that of an infinite circular cylinder 
filled with fluid. The boundary of the cylinder is a material shell which is assumed 
to rotate about its center of symmetry with angular velocity a — a(t) (positive 
rotation is counterclockwise). We restrict our attention to planar and circularly 
symmetric motions, so that the fluid velocity u is given by 

/ .s i ^ (-32,31,0) (-3: 2 ,2:i,0) 
U = U{X 1 ,X2,X 3 ,t) = V(r,t) = V(X 1 ,X2,t) , 

r r 

with r = \J x\ + x\. This symmetry assumption is consistent with the Navier- 
Stokes equations as long as the initial velocity is planar and circularly symmetric, 
i.e., weak solutions satisfying the symmetry assumptions exist globally in time for 
symmetric initial data. It should be remarked that the assumption that such flows 
stay planar for all Reynolds number is unphysical - 3D turbulence is expected if 
the Reynolds number is large enough. The limit behavior we wish to present is 
motivated by mathematical considerations. 

Let us first observe that, if velocity is circularly symmetric, and, in particular, 
tangent to concentric circles around the origin, vorticity is also circularly symmetric, 
and therefore constant in the same circles. This implies that Vw is perpendicular 
to these circles, and therefore u and Vw are orthogonal everywhere. Therefore, the 
nonlinearity u ■ Vcj is the vorticity equation vanishes identically, and the vorticity 
equation becomes the heat equation. This basic fact is what it makes possible to 
prove the following result. 

Theorem 4. Let u v be the solution of the 2D Navier-Stokes equations in the unit 
disk, with no slip boundary data with respect to boundary rotation with prescribed 
angular velocity a G BV(M). Assume that the initial velocity u G L 2 (T>) has 
circular symmetry, i.e. Uq(x) — vq(\x\)x . Then, uq is a steady solution of the 2D 
Euler equation and u v converges strongly in L°°([0, T], L 2 (T>)) to uq as v — ► 0. 

The proof of this result is based on a semigroup treatment of a symmetry-reduced 
problem. The important issue is that the symmetry assumption not only prevents 
boundary layer separation, but it also eliminates the nonlinearity of the problem. 
Thus, the symmetry-reduced problem is a linear, singular coefficient perturbation 
of the heat equation, and its treatment is classical. 
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Let D = {\x\ < 1} C M. 2 . We begin by considering the following 2x2 system of 
PDE, consisting of a pair of uncoupled heat equations: 



Let us suppose that uq is covariant under rotations, i.e., uo(Rgx) — RqUq(x) for 
any rotation Rg. Then (|4"3")l is the symmetry-reduced 2D Navier-Stokes system of 
equations with initial velocity uq- 

Let us assume, to begin with, that uq € C°° (D) is covariant under rotations and 
that uq and a <E C°°(K) satisfy the following compatibility conditions: 



Then it follows that there exists a unique weak solution u v — u v (x, t) e L 2 (R; H 4 (D))C) 



H 1 (R;H 2 {D)) n H 2 (R;L 2 (D)) of (03j In particular we have u v G C°(R; C l {D)). 



With this we can now establish the following lemma. 

Lemma 5. Let uj v = V • u v = d Xl u^ — d X2 u\. Then uj u is a classical solution of 



Additionally, J D oj u {x, t) dx — a(t) at every t > 0. 

Proof. Start by noting that, if uo is covariant under rotations about the origin, then 
luq is circularly symmetric and, hence, so is uj u (-,t); this is due to the rotational 
invariance of the heat equation. 

The partial differential equation in ([45]) can be trivially deduced from |43|) , along 
with the initial condition, so that all that remains is to verify that the boundary 
condition is correct. 

We integrate (|4"S"|) in space and use the divergence theorem to obtain: 



(43) 




(44) 



Uq{x) — a(0)§^r on 3D; 
Au =a'(0)f^ on 3D. 



(45) 




in (0, oo ) x D; 

on {\x\ = 1} x [0, oo) 

in D x {t = 0}. 





where we used the circular symmetry of io v in the last step. 
Next observe that lo v = — div (u") 1 - , so that 




as u v is a solution of (|43| . 

These two facts yield the desired conclusions. 



□ 
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Let us note in passing that the identity J D uj^(x,t) dx — a(t) is valid (a.c. in 
time) under the much weaker assumption uq G L 2 {D) and a G BV(M.), without the 
compatibility conditions (f4"4")) , since then the solution u v G L 2 ((0, T); 

Next we state and prove our main theorem in which we examine the inviscid 
limit of ui v . 

Theorem 6. Let a G BV(M) and assume that a is compactly supported in (0, oo). 
Let Uq G L 2 (D;M. 2 ) be covariant under rotations and assume that uiq = V ■ Uq G 
L 1 (D). Then we have: 

(1) There exists a constant C > such that 



Jd 

(3) We also have, for any < a < 1, 

\W - Wo||z,i({|x|<a}) -> as v -> 0, 

a.e. m iime. J/wo € C°(D) then the convergence is uniform in time. 

Given a G BV(M), we call the modified BV function 5 the left-continuous 
"correction" of a. 

Proof. We start by choosing appropriate regularizations of a and of uq. Let a n be 
the regularization obtained in Lemma 4.1 of [19] , satisfying the following properties: 

(1) a n G C°°(R); 

(2) a n — > a pointwise and strongly in Lj oc (R), for any 1 < p < oo; 

(3) a' n a' weak-* in BM(K). 

Above, the derivatives are understood in the sense of distributions. Let ito, n be a 
sequence of vector fields in C%°(D), commuting with rotations, such that Uo,„ — > uo 
strongly in L 2 (D). Since a is compactly supported in (0, oo) it follows immediately 
that the compatibility conditions (14411 are satisfied for ito, n and a n . Hence, if u v n 
denotes the solution of the symmetry-reduced Navier-Stokes equations (|43|) with 
data up,™ and a n then io v n = V • is a solution of (|45p such that 



by Lemma [5] 

Let us prove the estimate in item (1) for u>^. 

We will derive an energy- type estimate in L 1 after appropriately regularizing 
|cj^|(a;,t) = sgn(ujn(x, t))uj!^(x, t). We first mollify y i— > |y|. For each e > there 
exists a C , convex function </> e , such that \4>' e \ is bounded, uniformly in e, and 
e — * | • | as £ — ► 0. 

Next, we multiply the PDE in |45|) by 4>' e (uj%) and integrate over the disk to 
obtain: 




for almost all time. 
(2) For almost every t G M, 




(46) 
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MOdx < " / A[0 £ «)]dx = v \ V[<t> s «)]-*dS(x) 

\x\<l «'|a|<l J\x\=l 

= <(t) 

2tt 



Since (j> £ E C 1 and convex, $>" e {yi v n ) > so that, by the divergence theorem 
d_ 
dt 

/ <p' e (^ n )dS(x)<C\a' n (t)\, 

l\x\ = l 

where we used the boundary condition for tu^ in (|45[) . Then, integrating in time 
gives: 

4> £ (uj u n ) dx < \ \oL n {j)\dr + / (j> e (u>o,n)dz. 

I\x\<l JO «'|a;|<l 

Observe that we can choose <p e in such a way that <j) e {y) — > \y\ monotonically. 
Consequently, by the Monotone Convergence Theorem, passing to the limit e — > 
in the expression above gives, for every < t < oo: 

(47) / \u v n (x,t)\dx< [ \a' n (T)\dr+ [ \uo, n (x)\ dx < \\a\\ B v + \\oo\\v, 

J\x\<l JO J\x\<l 

from which we obtain item (1) for ui^. 

We have shown that, for each v > 0, {^n}^=i i s uniformly bounded in L°°(R; 
Thus, passing to subsequences as needed, we find that {uj^}^ =1 is weak-* L°° (R; BM(D j) 
convergent as n — > oo. 

On the other hand, the properties of a n , together with the proofs of Proposition 
4.3 and Proposition 5.1 in [TU] imply that u v n — u v weakly in L 2 ((0, T) x D), where 
u v is the solution of (|4"3")l . with data uq and a. It follows immediately that 

< = V x • U ", 

weakly in £ 2 ((0, T); H~ 1 (D)). Uniqueness of limits implies that the convergence of 
Lo v n to u> u can be improved to weak-* L°°(R; BM.(D)). We use again an estimate 
obtained in the proof of Proposition 5.1 of [T5], see (5.3) there, to find that {w^} n 
is bounded uniformly in L 2 ((0,T) x D). Hence we obtain 

uo v n w v weakly in L 2 ((0, T) x D). 

Therefore, since L 2 ((0,T) x D) C i 2 ((0, T); L 1 (£))), it follows that w"(-,t) G 
L 1 (Z>) a.e. t > 0. Weak lower-semicontinuity of the L 1 -norm implies that item (1) 
follows from (|4"7) . 

Since a„ — > a strongly in i 2 Qc (R) it follows that, for any <p = (f(x) G L°°(D) 
and ^ = ^(*) € i 2 ((0, T)) we have 

- t)^ fp(x)if)(t) dxdt — jT y f ^ - w" (x, i)^ dxd*. 

Take y(x) = 1. Then, using (HBJ), we see that the left-hand-side above vanishes 
identically, so that 

J (^Q-LJ V (x,t)j<lp(t)dxdt = 0, 

i.e. 



a(t)-J LJ v (x,t)dx\ ip(t)dt = 
for any ip G L 2 ((0, T)). Take V = ip(t) = a(t) - / D u v (x, t) dx; this gives item (2). 
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Note that cd" is a solution of the heat equation with viscosity v in D, with initial 
data ujq. 

To prove (3), let £1 be an open, compactly contained subset of D and let <p E 
C^°(D) be such that <f> = 1 in a neighborhood of Ct. We consider v v = <\><jJ v , 
extended to the full plane. Then v v satisfies: 

d t v v = uAv u +F u , 

with 

F v = -u(2V<f> ■ Vuj v + uj v A<j)). 
We apply Duhamel's formula in the whole plane to obtain: 

(48) v v =e WA (0w o )+ / e v( - t -^ A F v (s)ds. 

Jo 

We wish to estimate v" = v v (x,t) for x £ Q. 

Clearly, the first term converges to uq in f2, as v — > 0, in whichever space ujq lies 
in. What we are left to prove is that the other term in (j48|) converges to zero. This 
follows from the fact that the heat kernel is sharply localized when v — > 0, together 
with the fact that the support of F v is bounded away from O, uniformly in v, as 
we are only interested in x G f2. 

□ 

From this result we derive three conclusions. The first is that, if a S BV(M.) 
then oj v — W weak-* in ^^ c ((0, oo); BM(D)) and hence W = luq + fJ>, where 
fi is a measure supported on the boundary of the disk D. This implies that, in 
considering the classical open problem of the inviscid limit for the Navier-Stokes 
equations in domains with boundary, at the very least one has to deal analytically 
with regularity at the level of vortex sheets, without a priori sign conditions. The 
second observation is that the vorticity generated by a boundary layer appears to 
be proportional to the acceleration of the boundary with respect to the adjoining 
flow. Third, we have proved that, if a is not constant, then the vorticity u> v does 
not converge in L 1 to the vorticity of the limit flow u>o- This shows that the L 2 
convergence of velocity fields obtained in |19j cannot be improved to convergence 
in derivatives. Finally, we would like to mention that, in [20] . we weaken the 
hypothesis that a € BV to a € L p , for p > 1, and we adjust the convergence 
obtained accordingly. However, the vorticity accounting we have shown here only 
works for a G BV. 

8. Conclusion 

As a conclusion for these notes, it would be interesting to highlight a few main 
points of the our discussion on the mathematical theory of boundary layers. 

(1) We remark on the extent of the current lack of understanding of the van- 
ishing viscosity limit in the presence of boundaries. At this stage, the 
Boundary Layer Problem formulated as formulated in Section 4 is wide 
open. 

(2) The Prandtl equation and its refinements take a magnifying glass approach 
to the difficulty involved in boundary layers. This approach is very interest- 
ing and useful for specific problems, but the difficulties in the well-posedness 
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of Prandtl's equation do not make it very promising as an avenue for the 
Boundary Layer Problem. 

(3) The magnifying glass approach works well for the case of Navier friction 
condition mainly because the boundary layer equation in this case is linear 
and well-posed. 

(4) At best, the Boundary Layer Problem involves treating a sequence of ap- 
proximate solutions to the Euler equations with vortex-sheet level regular- 
ity, which is the critical case for the weak continuity of the nonlinearity. 

(5) The strength of the boundary layer as a vortex sheet, and therefore some 
of the difficulty in the Boundary Layer Problem appears to be associated 
with the acceleration of fluid past the boundary. 

Finally let us put together a few problems connected with what was presented 
here. First, establish rigorously the ill-posedness of Prandtl's equation without 
monotonicity. Second, we mention extending the LP control on the vorticity in the 
case of Navier friction condition to p < 2. Third, in [22], the author, together with 
H. Nussenzveig Lopes and Zhouping Xin introduced the notion of boundary coupled 
weak solution of the incompressible Euler equations, and we proved the existence 
of such a weak solution in the case of the half-plane. It would be interesting to 
know if the vanishing viscosity limit, for example, for Navier friction condition in 
the half-plane, gives rise to such weak solutions. Also, if vorticity is L p , p > 2 the 
theory of renormalized solutions of DiPerna-Lions implies that weak solutions of 
the Euler equations conserve p-norms exactly, for flows in domains without bound- 
ary. Does this remain true in the case of vanishing viscosity limits in domains with 
boundary? Would it be possible to find an example, in the spirit of the circularly 
symmetric flows in Section 7, but for which there is boundary layer separation? One 
possibility is to look for the solution of the Navier-Stokes system on the disk, with 
initial vorticity given by an odd eigenfunction of the Dirichlet Laplacian. Finally, 
is boundary layer separation possible for flows with Navier friction condition? A 
related problem is to study the vanishing viscosity limit with Navier friction con- 
dition in nonsmooth domains. Another pair of problems is to extend the boundary 
layer analysis done in Section 6 either to p = oo or p < 2. 
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